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ABSTRACT 

Grain growth by the accretion of metals in interstellar clouds (called 'grain growth') could 
be one of the dominant processes that determine the dust content in galaxies. The impor- 
tance of grain size distribution for the grain growth is demonstrated in this paper First, we 
derive an analytical formula that gives the grain size distribution after the grain growth in 
individual clouds for any initial grain size distribution. The time-scale of the grain growth 
is very sensitive to grain size distribution, since the grain growth is mainly regulated by the 
surface to volume ratio of grains. Next, we implement the results of grain growth into dust 
enrichment models of entire galactic system along with the grain formation and destruction 
in the interstellar medium, finding that the grain growth in clouds governs the dust content 
in nearby galaxies unless the grain size is strongly biased to sizes larger than ^ 0.1 |i.m or 
the power index of the grain size distribution is shallower than ^ —2.5. The grain growth 
in clouds contributes to the rapid increase of dus t-to-gas rati o at a certain metallicity level 
(called critical metallicity in Asano et al.ll201l1 and llnouel201 II) . which we find to be sensitive 
to grain size distribution. Thus, the grain growth efficiently increase the dust mass not only 
in nearby galaxies but also in high-redshift quasars, whose metallicities are larger than the 
critical value. Our recipe for the grain growth is applicable for any grain size distribution and 
easily implemented into any framework of dust enrichment in galaxies. 
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1 INTRODUCTION 

In the interstellar medium (ISM), dust grains are the most efficient 
absorber of stellar light. The spectral energy distributions and the 
radiative heating and cooling in gala xies are thus stron gly regulated 
by the presence of dust (e.g. Yamasawa et al.ll201 ih . This means 
that the understanding of dust enrichment in galaxies is crucial in 
the studies of galaxy evolution. 

Dust enrichment in galaxies is gover ned by vario us pro- 
cesses depending on age, metallicity, etc. jDwekl [l998h . In the 
earliest stage of galaxy evolution , dust is predomin antly pro- 
duced by supernovae (SNe) (e.g. iKozasa et alj |2009|) , while at 
later epochs asymptotic giant branch (AGB) stars also contribute 
jValia nte et al. 2009). The ti me-scale of dust destruction by SN 
shocks is a few x 10* vr J Jones. Tielens. & HoUenbachl [l996l : 



ISerra Dfaz-Cano & Jone sl l2008h . while that of dust supp ly from 
stella r sources is longer than 1 Gyr in the Milky Way dMcKed 
Il989l) . Therefore, dust grains should grow in the ISM by the accre- 
tion of metals onto grains (Draine 2009) to explain the existence of 
dust in the ISM. The growth occurs most efficiently in molecular 
clouds, where the typical number density of hydrogen molecules 
is ~ 10 cm" jHirashitallioOOah . This process is called 'grain 
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growth in clouds' in this paper. Observational pieces of evidence 
for the grain growth in clouds come from large r depletion of metal 
elements in cold clouds than in warm medium dSavage & SembachI 
1 19961) . 

A lot of chemical evolution models treat the evolution of dust 
content in galaxies. These models usually include dust produc- 
tion by stars, grain grow th in clouds, and d ust destruction by SNe 
(e.g. lDweklll998l: llnoud 120 03; Zhukovska . Gail. & Trielo3l2008l: 
ICalura. Pipino. & Matteuccili2008. : i Asano e~.ll201 ih . Most of the 
models which consider the grain growth in clouds indicate that this 
process dominates the dust budget at sub-solar or solar metallici- 
ties. The grain growth occurs through the accretion of metals, so 
that the increasing rate of grain mass by the accretion of metals is 
proportional not only to the metallicity but also to the grain surface- 
to- volume ratio, which is very sensitive to the grain size distribu- 
tion. 

Most models so far assume a certain grain size distribution or 
a typical grain size to estimate the surface-to-volume ratio. How- 
ever, since the dominant processes governing the grain size distri- 
bution shouMj^arywi^^ Mathi| 
1 19971 : iHirashita et alJbOld : lYamasawa et al.ll201 ih . it is expected 
that a variety of grain size distributions emerge in a complex way 
depending on age and metallicity. The first source of dust in the 
history of galaxy evolution is SNe, and the dust grains produced by 
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SNe are possibly biased to large (> 0.1 \im) sizes because small 
grains tend to be destroyed in the shocked r egion within SNe be- 
fore being injec ted into the interstellar space (Bianch i & Schneidej 
l2007l:lNozawae t al. 2007). Hirashita et al. ( 2010) show that small 
grains are produced by shattering driven by interstellar turbulence 
if the dust abundance i s as h igh as that expected from the solar 
metallicity. Ijones et al] ( Il996h show that shattering in interstellar 
SN shocks increases the abundance of small grains. Efficient pro- 
duction of small grains enhances the surface-to- volume ratio, acti- 
vating the grain growth by the accretion of metals. Thus, we should 
consider various grain size distributions depending on galaxy age 
and metallicity, and the grain growth efficiency may vary with a 
variety of grain size distributions. 

The first aim of this paper is to formulate the grain growth 
in clouds by explicitly considering the dependence on grain size 
distribution. Therefore, the former part of this paper is devoted to 
the formulation of the grain growth in clouds under an arbitrary 
grain size distribution. Then, by using this formulation, we point 
out the importance of grain size distribution for the grain growth by 
accretion. The final scope of this work is to examine if the grain size 
distribution has a significant influence on the dust enrichment in 
galaxies through the grain growth in clouds. Thus, in the latter part 
of this paper, we implement our formulation of the grain growth in 
clouds into a simple framework of dust enrichment in a galaxy, also 
taking into account the grain formation by stellar sources and the 
destruction by interstellar shocks driven by SN remnants. Thereby, 
we will show that our formulation of dust growth is successfully 
incorporated into dust enrichment models, and we will address the 
importance of grain size distribution for the grain mass budget in 
galaxies. 

This paper is organized as follows. We explain the formula- 
tion in Section (2] and describe some basic results on the evolu- 
tion of grain size distribution through the grain growth in individual 
clouds in Section[3] We implement the results for the grain growth 
in clouds into a simple evolution model of dust mass in an entire 
galactic system in Section |4] where the model also treats the dust 
formation by stellar sources and the dust destruction by SN shocks. 
We discuss the results in more general contexts in Section |5] Fi- 
nally, Section|6]gives the conclusion. 



2 FORMULATION 

In this section, we formulate the evolution of grain size distribution 
by the accretion of metals (grain growth) in a single interstellar 
cloud. Our procedures in this paper are divided into the following 
two steps: (i) we construct a formulation of the grain growth in 
clouds, which is conveniently incorporated in any dust evolution 
models in galaxies; and (ii) we show that our formulation in this 
section can be used generally to treat the grain growth in galaxies. 
This section is aimed at item (i). In Section|4l we address item (ii) 
by incorporating our formulation for the grain growth into simple 
dust enrichment models which also include dust formation by stel- 
lar sources and dust destruction in SN shocks. Other mechanisms 
that modify the grain size distrib ution such as shattering and coagu- 
lation are treated in other pape rs ( Jones et al.l 19941 19961: lYan et all 
l2004l : lHirashita&Yai]|2009l : lYamasawa et al.ll201lh~ These pro- 
cesses are to be included in future work for the comprehensive un- 
derstanding of the evolution of grain size distribution. 

Throughout this paper, we call the elements composing grains 
'metals'. We only treat grains refractory enough to survive after 
the dispersal of the cloud, and do not consider volatile grains such 



as water ice. More specifically, we consider silicate and graphite as 
main dust components. We also assume that the grains are spherical 
with a constant material density s, so that the grain mass m and the 
grain radius a are related as 
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(1) 



2.1 Evolution of grain size distribution 

We define the grain size distribution such that n(a, t) da is the 
number density of grains whose radii are between a and a + da 
at time t. For simplicity, we assume that the gas density is constant 
and the evolution of grain size distribution occurs only through the 
accretion of metals on dust grains. In this situation, the number den- 
sity of grains is conserved. Thus, the following continuity equation 
in terms of n(a, t) holds: 



dn(a, t) d , , , , 



0, 



(2) 



where u = da/dt is the growth rate of the grain radius and is 
given in the next subsection. Coagulation and shattering, which do 
not conserve the number density of grains, are not treated to focus 
on the grain growth by accretion here. Note that these processes 
do not change the grain mass, while the grain growth by accre- 
tion increases it. The evolution of grain size dis tribution by coag- 
ulation and shattering is treated in other papers jjones et alJI 19961 : 
iHirashita & Yaij2009l:lo"rmel et al.l2009h . We also neglect possible 
grain destruction mechanisms in molecular clouds by cosmic rays 
or shocks. 



2.2 Grain growth rate 

The grain growth rate is basically determined by the collision 
rate between a grain and particles of the relevant metal species. 
We adopt silicate and graphite as dominant grain species (e.g. 
iDraine &Le3ll984 . and denote elements composing the grains as 
X (for example, X = C, Si, etc.). We neglect the effect of Coulomb 
interaction on the cross section (i.e. the cross section of a grain 
for accretion of metals is simply estimated by the geometric one) 
because the grains and the a t oms are neutral in molecular clouds 
JWein gartner & Drainj|l999t lYan et _a l."2004'). In fact, the ioniza- 
tion degree in dense clouds is 10~® (|Yan et al. 2004), which means 
that almost all the metal atoms colliding with the dust grains are 
neutral. The rate at which atoms of element X strike the surface 
of a g rain with radius a is denoted as TZ and is estimated as jEvanj 
I1994I) 

n = 4naW(^^y^\ (3) 
V 27rmx / 

where nx is the number density of element X, k is the Boltzmann 
constant, Tgas is the gas temperature, and mx is the atom mass 
of element X. In general, dust grains are not composed of a sin- 
gle species. We adopt a key element, whose mass fraction in the 
grain material is /x, and represent the grain growth by the accre- 
tion of element X. The concept of key element is also adopted by 
IZhukovska et all feOOSi) . 

By using fx, the increase of the grain mass m is estimated as 

dm 1 

= /x mxSTZ, 



(4) 



where 5* is the sticking probability. This equation is converted into 
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the increasing rate of a by using equations l[T) and l|3j as 



da 



nxft) mxg / fcTgas \ 



1/2 



(5) 



/xs V27rmx/ 

In fact, nx (t) is a function of time because the metal abundance in 
gas phase decreases as the dust grains grow. This effect is treated 
in the next subsection. 



2.3 Depletion of gas-phase metals 

The decreasing rate of the number density of element X in gas 
phase is equal to the grain growth rate per volume: 



^ = - /" nSn{a, t)Aa. 
™ Jo 

Now we introduce the ^-th moment of a as 



(a )(t) = — / a n{a, t) da, 
rid 



(6) 



(7) 



where is the number density of dust grains, which is indepen- 
dent of t: 

poo 

rid = n{a, t) da. (8) 
Jo 

The moments are functions of t, and their values at i = are de- 
noted as (a*)o {= (aO(0))- By using the second moment of a and 
equation ([3} for TZ, equation (|6} can be expressed as 

Note that the third moment of a is related to the dust mass 
density pd{t) as 

Pd(i) = ^7r(a^>(t)snd. (10) 

The initial dust mass density is Pd(0) = |7r{a^)o sud- 

Here we quantify the initial number density of element X. The 
total number density of element X both in gas and dust phases is 
written as 



(9) 



(11) 



where Z is the metallicity, and (X/H)q is the solar abundance rel- 
ative to hydrogen in number, and nn is the number density of hy- 
drogen nuclei. We denote the initial fraction of element X in gas 
phase as ^: 

nx(0) = ^nx,tot. (12) 

Since the initial number density of element X in dust phase is 

/jfPd(0)/mx, 

/xPd(O) 



mx 



(1 — 5) nx,tot. 



By using equation dlOt at f = 0, equation ( I13t is written as 



1-C 



|7r(a^)o/xsnd 
mxnx.tot 



(13) 



(14) 



Thus, the normalization of the grain size distribution is determined 
by 



rid 



mx(l-e) Z\ (X 



|^(a3)o/xs \7^q) (h)© 
where we used equation dl It for nx,to 



n-H, 



(15) 



2.4 Formal solution 

The most important characteristics of the grain growth by accretion 
is that the increasing rate of a is independent of a (equation |5]l- 
Thus, equation ^ gives a formal solution as 



n(a, t) = n(a - A{t), 0), 
where 



~ /xs \27rmx/ 



nx{t) dt. 



(16) 



(17) 



A{t) is obtained if we give nx{t) by using equation (|9]l. We for- 
mally assume that n(a, 0) = for a < so that equation l |16l l can 
be used even for a — A{t) < 0. 

We introduce the following indicator for the dust mass in- 
crease in the cloud (equation [TO): 

PdW ^ 

Pd(0) {a3>o 

(a^)o + (a^)oA(t) + {a)oA{tf + Ajtf 



(18) 



where we expand the third moment by using equation ( |16t . Note 
that the dust mass in the cloud at t is {a^){t) / (a^)o times the initial 
dust mass. A similar equation is also obtained for the grain mantle 
growth as shown bv iGuillet. Pineau des Forets. & Jones, C2007i) . 



2.5 Typical time-scale 

For numerical calculations and interpretations of the results, intro- 
ducing a typical time-scale is convenient. The typical time-scale of 
grain growth by accretion is defined by 



do 



mxnx.totS ( kT. 



/xs \27rmx/ 



1/2- 



(19) 



where ao is a typical grain radius given arbitrarily. By using r, 
equation ( 117b is reduced to 



A{t)^ 



ao 



nx{t) dt, 



(20) 



(21) 



while equation l|9ll is written as 

dnx(t) _ 3nxao , _ (a^)ft) 
dt - r ^' (a3)o ' 

where we have used equation ( I14t . Note that 

{a^){t) = (a2>o + 2{a)oA{t) + A{tf. (22) 

A set of equations (|20M22t is solved to obtain the grain growth 
Ait). 

2.6 Selection of quantities 

We co nsider silicate and g raphite as representative grain compo- 
nents ( lDraine&Le3ll984l) . In order to avoid the complexity aris- 
ing from compound species, we treat those two species separately 
as a first approximation. The quantities adopted in this paper are 
sum marized in Tab l e [T] 

Ijones & NuthI ( 1201 ih point out that the accretion of Si, Fe, and 
Mg in the presence of abundant H2, CO and H2O may lead to for- 
mation of a complex ice. Therefore, the real picture of the silicate 
growth in clouds may be chemically complicated. They also state 
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Table 1. Adopted quantities. 



Species 


X 


/x 


mx [amu] 


(X/H)0 


s [g cm 3] 


Silicate 


Si 


0.166 


28.1 


3.55 X lO"-'^ 


3.3 


Grapliite 


C 


1 


12 


3.63 X 10-"' 


2.26 



that any silicate materials produced in the ISM with reasonable sce- 
narios do not show the spectral properties actually observed in the 
ISM. For carbonaceous dust, graphite and amorphous carbon have 
different chemical properties, leading to different growth properties 
in accretion. In this paper, because the knowledge about chemical 
properties is still poor, we simplify the picture by assuming that the 
grain growth is regulated by the sticking of the key species. 



2.6.1 Silicate 

We assume that the accretion of silicon regulates the growth of 
silicate grains. A larger abundance of oxygen is compensated by 
a lar ger number required to compose silicate. IZhukovska et al.l 
( I2OO81) also consider Mg as well as Si as a key species. The abun- 
dances of Mg and Si are almost the same if we assume the solar 
abundance pattern. Therefore, the following result does not change 
significantly even if we adopt O or Mg as a key element. 

The solar abundance of Si is (Si/H)Q = 4.07 x lO"-'^ 
( lLoddersll2003h . and the mass of a Si atom is msi — 28.1 amu 
(1 amu = 1.66 x 10"^^ g). For the material properties of silicate, 
we assume s = 3. 3 g cm~'^ and .fsi ~ 0.166 based on composition 



Mgi.iFeo.9Si04 jPraine & Leei[l984 . By using equation l[TT), the 
typical time-scale r (equation|I9t can be estimated as 

r = 6.30 X lO'ao.i(^/Z0)"'n3-^r3o''"S'o::; yr, (23) 

where ao.i = ao/0.1 |J.m, na = nn/lO^ cm^^, Tso = Tgas/SO K, 
and So. 3 = S/0.3. Unless otherwise stated, we adopt wh = 
10"^ cm~^ and T = 50 K for the typical values d erived from obser - 
vational properties of Galactic molecular clouds (Hirashita 2000a ), 
and S = 0.3 ( Leitch-Devlin & Williams 1985; Grassi et al. 20 lH). 
Although S may be almost I in such a low tempe rature environ- 
ments as molecular clouds JZhukovska et al.|[2008[) . the chemical 
factor (i.e. if the sticked atom finally becomes the p art of the grain 
or no t) is quite difficult to quantify as noted above jjones &Nutlj 
Thus, we adopt a conservative value of S in this paper. In 
fact, T depends on ZS, so a larger value of S is compensated by 
a smaller Z. In other words, the same dust growth rate is achieved 
with metallicity ^/(S/0.3)ifwe adopt a different value of S. Sim- 
ilar discussion also holds for the other quantities (T and n). 



2.6.2 Graphite (carbonaceous grains) 

The s olar abundance of C is (C/H)q = 2.88 x 10"* fcoddersl 
l2003h and the mass of a C a tom is mc = 12 am u. For graphite, we 
assume s = 2.26 g cm~^ toraine & Le3ll984l) and /c = 1 (i-e. 
graphite is solely composed of C). The typical time-scale is 

r = 5.59 X IQ'ao.iiZ/T.QY^n-^^T-^'-^^So^ yr. (24) 

If we adopt other carbon species such as amorphous carbon, the 
time-scale does not change as long as s is similar. If a light material 
like hydrogenated am orphous carbon (a-C:H) (s — 1.3 -1.5 g cm""^ 
and /c = 0.91-0.98: [Serra Diaz-Cano & Jonesll200i) is adopted, 
r becomes 0.52-0.65 times the above value. The face value of r 



indicates that a-C:H grains grow more efficiently than graphite, al- 
though the difference in S between these two species is not clear. 
The destruc tion of a-C:H is, how ever, also more efficient than that 
of graphite JSerra Dfaz-Cano & Jones..2008i) . 



2.7 Initial grain size distribution 

We examine a variety of initial grain size distributions. The various 
models surveyed are called Models A-H as listed in Table|2] 



2.7.1 Single size ( 5 function ) 

We first examine the case where the grain size distribution is 
strongly peaked at a certain typical radius. For example, dust 
grains condensed in SNe have a typical ra dius of ~ 0.1 \im if 
the shock efficiently destro ys smaller grains dsianchi & SchneideJ 
I2OO7I : iNozawa et al.ll2007l) . Here, the typical radius is assumed to 
be ao, and S function is adopted for the grain size distribution for 
simplicity. Then we obtain 



n(a, 0) = ndS{a — ao), 

where rid is given by equation il5i . The moments are (a*)o 



(25) 



ao. 



2.7.2 Power law 



If the grain size distribution is described by a power law, the typical 
grain size is not obvious. Thus, we arbitrarily adopt ao = 0.1 i^m 
for power-law grain size distributions. The upper and lower bounds 
for the grain radii are denoted as amin and Omax, respectively: 



n{a, 0) 



(1 - r)nd 

l-r _ „l-rN 
max "min/ 



ln(aniax/amin) 



{a- 



(r = l); 



for dm in ^ a ^ If a < amin or a > amax 

The normalization rid is given by equation (Tsj. 
The moments are given by 



(26) 

n{a, 0) = 0. 



(a*)o 



£ + 1- 

1 ^ 



Injg/p) 



-e 

amax 



(r/1, ^ + 1); 

tin 

- (r = l); 

ln(amax/amin) (r = £ + 1). 



(27) 



3 GRAIN GROWTH IN INDIVIDUAL CLOUDS 

First, we examine single size cases (Models A and B), where the 
grain size distribution is described by 5 function. For the initial 
fraction of element X in gas phase, we adopt ^ — 0.5 and 0.9 
as representative cases where a moderate fraction and only a small 
fraction of metals are in dust phase, respectively. In Fig.[T] we show 
the evolution of nx(t)/nx, tot and {a"^)(t)/(a^)o for Models A and 
B. Note that (a^)(t)/{a^)o is proportional to the dust mass or the 
dust mass density (equation 11 8t and that nx{t — 0)/nx,tot ~ 
^ (the fraction of metals in gas phase at t — 0). As the model 
imposes, the dust mass increases with the depletion of gas-phase 
metals onto the grains. For larger ^, the grains continue to grow to 
larger radii and for a longer time, since the available metals in gas 
phase is more abundant and the abundance of dust relative to the 
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Table 2. Models. 



Model 




ao 




r 


'^min 


Q-max 


{a)o 






(a3)o/(a2)o- 






[lim] 






[|xm] 


[|xm] 








[l^m] 


A 


5 


b 


0.5 
















B 


5 


b 


0.9 
















C 


P 


0.1 


0.5 


3.5 


0.001 


0.25 


0.00167 


0.00216 


0.00402 


0.0159 


D 


P 


0.1 


0.5 


2.5 


0.001 


0.25 


0.00281 


0.00667 


0.0158 


0.0887 


E 


P 


0.1 


0.5 


4.5 


0.001 


0.25 


0.00140 


0.00153 


0.00187 


0.00279 


F 


P 


0.1 


0.9 


3.5 


0.001 


0.25 


0.00167 


0.00216 


0.00420 


0.0159 


G 


P 


0.1 


0.5 


3.5 


0.0003 


0.25 


0.000499 


0.000659 


0.00156 


0.00874 


H 


P 


0.1 


0.9 


3.5 


0.003 


0.25 


0.00499 


0.00633 


0.0103 


0.0273 



Initial grain size distribution: "p" for the power law and "(5" for the 5 function. 
^ The time-scale is scalable for any oq by using equations j23t and i24\ . 
The moments of a are not free parameters but are calculated once the other parameters are fixed. 
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Figure 1. Evolutions of the normalized grain volume (or mass) relative 
to the initial value, (a'^) / (a^)(), and the number density of element X in 
gas phase normalized to the total (gas phase + dust phase) number density, 
^x/'ix tot. along with the normalized time t/r. The 5 function grain size 
distribution is adopted. Solid and dashed lines show the results of Models 
A and B, respectively (i.e. ^ = 0.5 and 0.9, respectively) with the upper 
and lower lines show the evolution of (a^) / (a'')o and nx/"x,tot : respec- 
tively. 



gas phase metals is smaller. At t S> t, {a^){t)/ {a'^)o approaches 
1/(1 — which means that the grain mass becomes 1/(1 — ^) 
times the initial mass after all the gas phase X is depleted onto the 
dust grains. 

Next, we examine the power-law grain size distri b utions 
(Models C-H in Table O. Mathi s. Rumpl. & NordsiecS h97l 
hereafter MRN) show that the extinction curve in the Milky Way 
can be fitted with a power-law grain size distribution. We fix the 
range of the grain size by adopting amin ~ 0.001 [im, and Omax ~ 
0.25 [im (MRN). We arbitrarily fix the value of ao as ao = 0.1 [im 
(Section r2.7.2t . Since the lower boun d of the grain size is poorly de - 
termined from the extinction curve dWeingartner & DrainelbOoH) . 
we also examine the case where the smallest radius is 0.0003 |J.m 
(3 A). We change the power-law index r from 2.5 to 4.5 to examine 
the cases where large and small grains are relatively dominated, re- 
spectively (note that r — 3.5 corresponds to MRN). We also vary 

In Fig. |2l we show the results for Models C-E. As shown in 
equation illl . the consumption rate of metals in gas phase is pro- 



Figure 2. Same as Fig.[T]but for the power-law initial grain size distributions 
(Models C-E). Solid, dotted, and dashed lines show the results of Models 
C, D, and E (r = 3.5, 2.5, and 4.5), respectively with the upper and lower 
lines show the evolution of (a^)/ (ii^)o ™d "x/'^x,tot> respectively. 



portional to the surface-to-volume ratio of dust grains. If r is large, 
the surface-to- volume ratio is large, so that the metals in gas phase 
is rapidly consumed onto dust grains. 

Fig. [3] illustrates the evolution of grain size distribution. To 
show the mass distribution per logarithmic size, we multiply a^ 
to n. Also in order to make the quantity dimensionless, we multi- 
ply ao/ud- We indeed observe that the small grains grow to larger 
grains. Since the increasing rate of grain radius is independent of a, 
the impact of grain growth is s ignificant for small grains as already 
shown by iGuillet et al.l ( |2007|) for the mantle growth. Moreover, 
gas-phase metals accrete selectively onto small grains (especially 
for larger r) because the grain surface is dominated by small grains; 
this is consistent with the results of Weinaartner & Draine ( 1999). 

To show the dependence on ^, we compare Models C and F. 
The results are shown in Fig.|4] Different values of ^ lead to differ- 
ent final-to-initial dust mass ratios (a'')/(a^)o — >■ 1/(1 — 1^), which 
are interpreted in the same way as in the 5 function cases (Fig.[T). 

Finally, in Fig. |5] we compare Models C, G, and H to exam- 
ine the dependence on amin- As expected, the results are sensitive 
to ttmin since the s urface-to-volume ratio o f grains differs largely 
(Table [2j (see also IWeingartner & Drainel | 19991). Therefore, it is 
important to specify the physics governing the smallest grain size 
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Figure 3. Evolution of grain size distribution for Models C, D, and E (Pan- 
els a, b, and c, respectively). The normalized times are shown in each panel. 
The size distribution is shown by ao(a/ao)*"(<ii *)/"d t" indicate the 
normalized grain mass distribution per logarithmic size. 



in interstellar clouds. Coagulation possibly depletes grains smaller 
tlian several A liHirashita & Yan.2009. : .Guillet et al..2011.) . 



4 IMPLEMENTATION INTO GALAXY EVOLUTION 
MODELS 

In this section, we relate tlie results above for the grain growth in 
individual clouds to the dust mass evolution in an entire galactic 



Figure 4. Same as Fig.|2]but for Models C and F = 0.5 and 0.9, respec- 
tively, with r = 3.5). 
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Figure 5. Same as Fig. |2]but for Models C, G, and H, where the minimum 
grain radius is varied as amin = 0.3, 1, and 3 nm, respectively, with r = 
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system. We adopt a chemical evolution model that calculates the 
metal and dust enrichment. To concentrate on the importance of 
grain size distribution, we simplify the galaxy evolution models, 
although our recipe for the grain growth is simple enough to be eas- 
ily implemented to more complicated galaxy evolution models. The 
entire galactic system is treated as one zone. The treatment of one 
zone is good if the dust, metals, and gas are mixed instantaneously 
or if we focus on a certain well mixed region in the galaxy. For 
example, in spiral galaxies which generally have a radial metallic- 
ity gradient, our models can be applicable to a certain radius range 
where the metallicity can be regarded as uniform. 

Since we are interested in the effects on dust enrichment, we 
compare our results with the dust abundance in galaxies. Given 
that the dust enrichment is closely related with the metal en- 
richment, it is convenient to derive the relation betwee n dust-to- 
gas ratio and metallicity (e.g. iLisenfeld & Ferraralll998l) . The dust 
mass is usually estimated by the far-infrared emission to trace 
the emission from large grains, which occupy a significant frac- 
tion of the dust mass. The dust mass estimated from the far- 
infrared emission can miss very small grains and very cold dust. 
The contri bution from very small grains to the dust mass is not 
significant jPesert. Boulanger. & Pugetll990l:lGalliano et al.l2005l: 
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ICompiegne et alj |201 ih . Very cold dust traced in longer wave- 
lengths than submillimetre may have a significant contribution to 
the t otal dust mass especially in metal-poor galaxies dOalliano et alJ 
12003 . 2005; Galametz et al. 2011), but its abundance is signifi- 
cantly affected by the assumed emissivity index of large grains. 
If we miss the contribution from very cold dust, the observational 
dust-to-gas ratio is underestimated in this paper, which enhances 
the importance of the dust growth in clouds to explain the addi- 
tional contribution from very cold dust. 

Although there is a large variety in the grain size distribu tion 
derived from the dust emission spectrum dOalliano et al.llioO^ . the 
regulation mechanism of grain size distribution in the ISM is not 
fully understood. Therefore, in this paper, we examine various grain 
size distributions (i.e. 5 function with various typical sizes and 
power law wit h various power i ndeces). Sputtering and shattering 
in SN shocks /jones et al.' '1996^) and shattering in interstellar tur- 
bulence ( Yan, Lazarian, & Draine 2004; Hirashita et al. 2ol3) 

may 

play a significant role in determining the shape of grain size dis- 
tribution. Inclusion of these processes into the evolution of grain 
size distribution is left for future work to concentrate on the grain 
growth in clouds in this paper. A simultaneous treatment of dust 
formation and destruction can be seen in Yamasawa et al. (201 1). 

Photo-processes are neglected in this paper. Dust may contain 
ice mantle when it is injected into the ISM from clouds. When the 
dust is exposed to stellar radiation, the evaporation of the ice man- 
tle may lead to grain disaggregation. If small refractory grains are 
ejected in the disaggregation, the total amount of refractory dust 
(silicate and graphite in this paper) does not change, so that our 
conclusion below is not affected. If the dust is destroyed and the 
refractory elements are returned into the gas phase in the photo- 
processing, the dust abundance may be affected by this process. 
The destruction of dust by photo-processing, if any, can be effec- 
tively included into the dust destruction efficiency (/?sn defined 
later), but we will see later that only the destruction by SNe is 
enough to explain the dust-to-metal ratio (Section|5]2j. 



4.1 Grain growth time-scale 

We define the increased fraction of dust mass in a cloud, /3, which 
is evaluated as 



A(Tci)[3(a')o + 3(a)oA(rci) + A(ta 



(28) 



where is the lifetime of clouds hosting the grain growth and 
equation l ll8b is used from the second to the third step. The dust 
mass in the cloud becomes (/3 + 1) times the initial value after 
the cloud lifetime, when the dust grown in the cloud returns in the 
diffuse ISM (see eguationfTSll. 

Now we consider how to implement our results into dust en- 
richment models for an entire galactic system. By using /?, the in- 
creasing rate of dust mass by accretion in clouds should be written 
as 



dMdust 



(29) 



where Mdust is the total dust mass in the galaxy, and is the 
mass fraction of clouds hosting the grain growth to the total gas 
mass. In equation l |29l ), we have assumed that the time-scale of dust 
enrichment is much longer than the cloud lifetime so that the mass 
increasing rate in each cloud is estimated by /3/rci times the dust 



mass in the cloud. On the other hand, iHirashital feOOOal) express 
the increasing rate of the dust mass in a galaxy through the grain 
growth in clouds by 



At 



(30) 



where Tgrow is the growth time-scale of the dust in a cloud, and ^ is 
the fraction of metals in gas phase (Section [2.2t . A similar expres- 
sion for the increasing rate of the dust mass is widely adopted in 
dust e nrichment models (e.g. lDwekll998l ; linouell2003l ; ICalura et al.l 
l2008h . 

Comparing equations ([29} and OOt . we obtain 



^grow ^ 



A(rei)[3(a2)o + 3(a)o^(rci) + A{t,,Y\ ' 



(31) 



where equation J28t is used for /3. Thus, we obtain Tgrow /tci as a 
function of t^\\ that is, if we give the lifetime of molecular cloud, 
we can evaluate Tgrow. Equation l |20t indicates that A(t^\) depends 
on T. Since t depends on the metallicity (equations 1231 and I24t. 
Tgrow /tci evolves as the system is enriched with metals. Here we 
adopt equation ([23} for t. To simplify the discussion, we focus on 
the metallicity and the grain size distribution with the other quanti- 
ties fixed; that is, t = 63.0ao.i (.Z/Zq)^^ Myr. As a result, if we 
specify Td and give a grain size distribution, we obtain Tgrow /tci 
as a function of Zr^. Thus, we can regard Tgrow /tci as a function 
of metallicity if we fix Td. The results are shown in Fig.|6] 

Fig- m shows that Tgrow /tci depends on the grain size distri- 
bution as expected from the results in the previous section. In par- 
ticular, if ao is smaller for the & function cases or r is larger in 
the power-law cases, the grains grow on a shorter time-scale so 
Tgrow becomes smaller. We also observe the dependence on if ^ 
is larger, the grains grow to a larger extent because the abundance of 
available metals relative to dust is larger. Thus, Tgrow /tci becomes 
smaller for larger ^ in the region of high Z or large Td. 

Fig- m also shows that the behaviour at small metallicities is 
described by Tgrow oc Z^^ . This is explained as follows. If the 
metallicity is low enough, Td ^ t. Thus, the decrease of gas-phase 
metals is negligible, and equation ( I20t can be written as follows if 
we recall that nx(0)/nx,tot = i'- 



A(Td) 2; C"OTd/T forTd<T. 



(32) 



By taking up to the first order for Tci/t, equation ( |31t is approxi- 
mated as 



^"'^^ -TOcZ-^ f0rTd«T. 



3(a2)o ao 



(33) 



This indicates two important properties. One is that the grain 
growth time-scale (which is proportional to Z~^) is roughly deter- 
mined by the typical growth time-scale t and the surface-to- volume 
ratio of dust grains. The other is that the cloud lifetime does not en- 
ter since the grain growth in each cloud during the cloud lifetime is 
slight enough. 

As Z becomes large, Tgrow /tci approaches a constant value 
independent of metallicity. Indeed, the grain growth saturates if 
the metals in gas phase are used out. Thus, (a'^)(Td)/(a^)o — >■ 
1/(1 ~ ^) as Z becomes large. In reality, the grain abundance in 
clouds would approach the value expected from the equilibrium be- 
tween the grain growth and the destruction by cosmic ray, shocks, 
ultraviolet light in clouds, etc. We expect that the destruction within 
the clouds does not affect the dust abundance in the entire galaxy 
since th e dominant de struction occurs in the diffuse ISM by SN 
shocks dMcKed [19891) . Thus, in the clouds, we only consider the 
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Figure 6. Ratio of the grain growth time-scale to the lifetime of cloud.s 
hosting the grain growth by accretion (rgrow/Td) as a function of Zt^\. 
(a) The cases of single grain sizes (<5 function size distributions) are shown. 
The solid and dashed lines show the cases for g = 0.5 and 0.9, respec- 
tively. The upper and lower Unes show the cases for ag = 0.1 and 0.01 
Jim, respectively, (b) The results for power-law grain size distributions are 
presented. The thick dotted, solid, and dashed lines are for r = 2.5, 3.5, 
and 4.5 with § = 0.5, while the dot-dashed line shows the case for r = 3.5 
with 5 = 0.9. The thin dotted lines show the relation predicted from equa- 
tion (33) (upper line) and equation | |3U with equation | |32) (lower line) (the 
dot-dashed line is relevant for comparison). 



grain growth and neglect dust destruction. In the dust enrichment 
model for the entire galaxy in Section |4!4] we include the dust de- 
struction by SN shocks in diffuse medium. By using equation i3U . 
we obtain 



(I-C)t-cI (forrci>r), 



(34) 



where we note that /3 ~ ^/(l — ^) for Td 3> t. This result indicates 
that the grain growth is regulated by the cloud lifetime since the 
dust grains use up all the gas-phase metals within Tci. The factor 
(1 — 5) indicates that the dust mass increases with a larger fraction 
if a larger part of metals are in gas phase. 

In Fig.|6j3, we overlay the relation predicted by equation (133 1 
for r = 3.5 (the upper thin dotted line, which should be compared 
with the dot-dashed line for the exact solution). Although it gives a 
good approximation for low metallicities, the discrepancy becomes 
significant at (Z/ZQ)(rci/10 Myr) > 0.1. This is because the 
increase of dust surface area by accretion further accelerates the 
accretion of gas phase metals. Such an effect of acceleration of 



grain growth is not taken into account in deriving equation ( 133b . 
but is included in the second and the third order terms of A in 
equation ( |31b . where A can be approximated by equation ( |32t if 
the saturation of grain growth by the depletion of gas-phase met- 
als is neglected. Thus, we expect that the approximation of Tgrow 
becomes better than equation ( 133) if we combine equations ( |31| l 
and ( I32t . The result predicted by this combination is also shown in 
Fig.|6]5 (thin dotted line). It fits the exact solution better than equa- 
tion ( 1331 ). but it underestimates Tgrow for large metallicities. This 
underestimate comes from the fact that we neglected the depletion 
of metals by the grain growth. This depletion effect occurs if the 
cloud lifetime is comparable or larger than the typical grain growth 
time-scale. Thus, we propose the following approximate formula to 
include the depletion effect: 



(35) 



where y = ao£,Tci/T (Equation [32}. In Fig.[7l we show this ap- 
proximate formula in comparison with the exact results. We ob- 
serve that the approximation is fairly good. By using equation J31| (. 
/3 — ^(rgrow/Tci)^^. Thus, we adopt the following approximation 
derived from equation ( I35t : 



1-C 



(36) 



with y = ao^Tci/r. This approximate formula has a merit that we 
do not need to solve the differential equation for the depletion of 
gas-phase metals (equation [9l(. Therefore, we hereafter use equa- 
tions i35l and \36i to estimate Tgrow and /3, respectively. 



4.2 Relation to the star formation 

Galaxies are enriched by metals and dust as a result of star forma- 
tion activities. Thus, it is usually convenient to model the quanti- 
ties in terms of the star formation. In fact, dust growth by accre- 
tion should be strongly related to star form ation, since b oth occur 
predominantly in dense clouds in galaxies ( lDweklll998h . Here we 
make use of /? and connect the star formation rate with the grain 
growth rate by accretion. 

If the fraction e of the cloud mass is converted into stars after 
the cloud lifetime Tci, the star formation rate, tp, is written as 



Tcl 



(37) 



where Afgas is the total gas mass in the galaxy. Then, equation i29\ 
can be written by using the star formation rate as 



dA'/, 



dust 



dt 

where D ; 



A/dust /A^gas is the dust-to-gas ratio. 



(38) 



4.3 Recipe for implementation of grain growth 

Here we summarize the recipe to include the grain growth by ac- 
cretion into dust enrichment models. There are two methods as de- 
scribed in sections 14.11 and 14.21 called Method I and Method II, 
respectively. 

Method I Assume Xd (the mass fraction of clouds hosting grain 
growth) and Td (the lifetime of these clouds). Calculate also (a*) 
for I — 1, 2, and 3 for a given grain size distribution (or use 
the values in Table |2}. Then, under any metallicity, r (the typical 
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Figure 7. Same as Fig.|6]but we show the approximate results with equation 
)35t . The thick lines show the approximate relations and, while the thin lines 
are calculated by the exact formulation used to produce the results in Fig.|6] 



grain growth time-scale) is calculated by equation ( 1231 for silicate, 
by equation J24b for graphite or by equation il9\ in general, and 
y = ao^Tci /r is also obtained. Note that ^ (the fraction of metals 
in gas phase) should be calculated within the framework of dust 
enrichment (Section l4.4t . Finally use equation OH (more conve- 
niently equation [35] for an approximate formula) to obtain Tgrow, 
which should be used in equation l |30t to obtain [dA/dust/dt]acc- 

Method II Assume Tci (the lifetime of clouds) and e (the star for- 
mation efficiency in these clouds). Calculate also (a*) for ^ = 1, 2, 
and 3 for a given grain size distribution (or use the values in Table 
O. Then, under any metallicity, r is calculated by equation i23l for 
silicate or equation ( I24t for graphite, and y = ao^Td/r is also ob- 
tained. Note that ^ should be calculated by a chemical enrichment 
model. Finally use equation l |28t (or more conveniently equation[36l 
for an approximate formula) to obtain /3, which should be used in 
equation l |38t . 



4.4 Dust enrichment model 

We demonstrate that our results so far are really implemented into 
galaxy evolution models with dust enrichment. The aim here is 
not to construct an elaborate chemical evolution model but to fo- 
cus on the effect of grain growth. W e adopt a one-zon e closed- 
box model of dust/metal enrichment bv lHirashital ( l2000bl) (see also 
lDwekll998l ; lLisenfeld & Ferrarall998[) . The model treats the evolu- 



tions of total gas, metals, and dust masses (Afgas, Mz, and Afdust, 
respectively) in the galaxy. In this model, the metals include not 
only gas phase elements but also dust. The equations are written as 



dt 

dMz 
dt 

dMdust 



dt 



-Zip + Ez 
-Vij + fi,,Ez 



+ 



dAfdus 



dt 



(39) 
(40) 
(41) 



where E and Ez are the rate of the total injection of mass (gas -1- 
dust) and metal mass from stars, respectively, /in is the dust con- 
densation efficiency of the metals in the stellar ejecta, and tsn is 
the time-scale of dust destruction by SN shocks. For simplicity, we 
do not treat individual element X but treat the entire metals; how- 
ever, if we are interested in a specific dust-composing element X, 
we can replace the relevant quantities with ones specific for ele- 
ment X. The simple treatments above are sufficient to demonstrate 
that our scheme of grain growth in clouds for any grain size distri- 
bution is really appl icable to the galaxy evolution models. We refer 
to othe r papers (e.g.! Dweklll998LIZhukovska et al.ll2008l:lGall et al.l 
l2011ah for more detailed dust enrichment models. 

Since we are not interested in the detailed history of dust 
production in stellar ejecta, we adopt the instantaneous recycling 
approximation; that is, a star with m > rrit (m is the zero-age 
stellar mass, and mt is the turn-off mass at age t) dies instanta- 
neously after its birth, leaving a remnant of mass Wm- Once the 
initial mass function (IMF) is fixed, the returned fraction of the 
mass from formed stars, TZ, and the mass fraction of metals that 
is newly produced and ejected by stars, Vz, are evaluated. Using 
these quantities, we write 



E = TZtp, 
Ez = {TZZ + yz)ip. 



(42) 
(43) 



We adopt 7^ = 0.18 and yx = 0.013 (Appendix). For /in, we 
examine two cases: /in = 0.1 and 0.0 1, which corr espond to the 
fiducial and the lower efficiency cases in^In oud ilOl ih . respectively. 

For the time-scale of dust destruction by SNe, we adopt an ex- 
pression rsN ~ A/gas/(esAfs7), where Es and Ms are the dust 
destruction efficiency and the gas mass swept by a single high - 
velocity SN blast, respectively, and 7 is the SN rate jMcKedl 19891) . 
Since we are interested in objects whose time-scale of star forma- 
tion is much longer than 10^ yr, it is assumed that the SN rate is 
proportional to the star formatio n rate (equation IA61 . We adopt 
esA/s = 1300 Mq jMcKedl 19891) . Then we obtain 



Me, 



TSN ~ 



(44) 



where /3sn = esMsj/ip ~ 9.65. As pointed out by I Jones & NuthI 
( l201lh . e is uncertain by a factor of 2 and is dependent on 
the as sumed grain composition (see also lSerra Dfaz-Cano & JonesI 
l2008h . 

Equations (|39M41I> are converted to the time evolution of 
the metallicity Z — Afz/Afgas and the dust-to-gas ratio 2? = 

Afdust /A/gas as 

Afgas dZ ^, 

d^=^^' 



i, 

Mgas dV 

1/) dt 



(45) 



- fin{nz + yz) - 

"dAfdust 



+ 



dt 



(46) 
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where we should evaluate [dMdust/dt]acc according to Method I 
or II in Section l431 It is convenient to combine the above two equa- 
tions to obtain the relation between V and Z: 



(/3sN + TZ)V 



1 



dMd 



dt 



(47) 



In Method I, [dMdust/df]acc/'(/' = I'C(f"SF/T'grow) (equa- 
tionl30t. where tsf = XdMgas/ip is the star formation time-scale. 
In Method II, [dAfdust/dtJacc/V' = /S^'/e (equation [38ll. A large 
rsF in Method I is equivalent with a small e in Methods II; that is, 
a small star formation efficiency means a long star formation time- 
scale. Because of this simple equivalence, we hereafter concentrate 
on Method II. Adopting Method II (i.e. equation (38]), equation l l47t 
is reduced to 

Vz^ = MTZZ + y z) - iPsN + 71- 13/ e)V. (48) 
aZ 



iLada. Lombardi. & AlvesI bOIO l) show that the star forma- 
tion efficiency in molecular clouds is roughly 0.1. They also men- 
tion that molecular clouds survive after the star formation activ- 
ity over the last 2 Myr. The comparison with the age of stellar 
clusters associated with mole cular clouds indicates that the life- 
time of cl ouds is ~ 10 Myr i Leisawitz. Bash. & Thaddeuslll989l : 
iFukui fel Cawamura 2010). Thus, we hereafter assume e = 0.1 and 
Tmoi ~ 10 Myr as standard values. For the initial condition, we 
assume D = and Z — 0. 

The relations between dust-to-gas ratio and metallicity are 
shown in Fig. [8] for fin = 0.1 and in Fig. |9] for fin = 0.01. We 
assume Zq — 0.015 l lLoddersl|2003l) . At low metallicity levels, 
the solution of equation \4Sl is approximated by I? ~ finZ. This 
is why the dust-to-gas ratio in the low-metallicity regime is lower 
in Fig. |9] than in Fig. [8] Above a certain metallicity level, a rapid 
increase of dust-to-gas ratio occurs because of the grain growth in 
clouds. The metallicity level at which this growth occurs is very 
sensitive to the grain size distribution. 

The observational data of nearby galaxies are also shown 
for comparison. For the uniformity of data, we select the 
samples o bserved by AKARl: blue compact dwarf galaxies in 
Hirashita & Ichikaw i ( l2009l) and spir al galaxies (M 8 1 from 



Sun &Hirashitd 1201 ll and M 101 from ISuzuki et al.l l2007t) . The 



dust masses are estimated from 90 \m\ and 140 yaa data by us - 
ing the mass absorption coefficient in lHirashita & Ichikawal ( l2009b . 
For the gas mass of the dwarf galaxies and M81, we adopt the 
H I mass, since the molecular mass is negligible or not detected. 
The H I masses of M 81 an d of the dwarf galaxie s are ta ken from 
IWalter et all ( |2008|) and in iHirashita & Ichikawal ( |2009|) . respec- 
tively. For M 101, we a dopt the sum of ne utral and molecular 
gas masses compiled in ISuzuki et al.l | |2007|) . The oxygen abun- 
dance is adopted for the indicator of the metallicity, and the solar 
abundance is assumed to be 12 + log(0/H) = 8.69 jLoddersI 
l2003h . The oxygen abundances of the dwarf galaxies are com- 
piled in|Hirashita & Ichikawa (2009), and t hose of t he spi ral galax- 
ies at the half-light radius are taken from iGamettI (l2002k . Typical 
errors of the observational quantities are comparable to the size 
of symbols in the figures. To ensure that the AKARI results are 
not systematically different from other results, we add the data 
in|lssa, MacLaren, & Wolfendale ( 1990), who estimated the dust 
mass from the extinction. The overall trend of the data is repro- 
duce d by the models. 

lAsano et al.l ( 1201 ih also include a sample whose dust mass 
is derived from Spitzer observations bv lEngelbracht et al.l ( l2008h 
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Figure 8. Relation between dust-to-gas ratio 23 and metallicity Z for var- 
ious grain size distributions with /;„ = 0.1. (a) 5 function grain size 
distributions (the solid and dotted lines present the results for ao = 0.1 
and 0.01 i^m, respectively), (b) Power-law grain size distributions (the 
solid, dotted, and dashed lines show the results for r = 2.5, 3.5, and 
4.5, respectively). The filled squares and asterisks represent the observa- 
tional data for dwarf and spi r al galaxies observed by AKARI, respectively 
fairashita & Ichikawd |2009| ; |Suzuki et al' 2007"; Sun & Hirashitd l201lh . 
The open diamonds are taken from Issa et al. ( 1990) as a spiral galaxy sam- 
ple. The cross on each line marks the critical metallicity. 



and show that the relation between dust-to-gas ratio and metallicity 
does not significantly change. The inclusion of Herschel data may 
boost the dust abundance because of the possib le contribution from 
very cold dus t especially in dwarf galaxies (i Grossi etal.ll2010h . 
iGalliano et al. I (^005) also find a large contribution from very cold 
dust in the submillimetre for some dwarf galaxies. However, since 
the modeling of submillimetre emission may be significantly af- 
fected by the assumed emissivity index of large grains, we do not 
use the submillimetre data in this paper. We should keep in mind 
that that inclusion of submillimetre data can rather decrease the 
dust mass especially for dust-i ich galaxies because th e dust temper- 
ature estimate becomes better dGalametz et al.l201 ll) . In their mod- 
els, the Spitzer 160 \im data are used; thus, it is not still unclear if 
there is a discrepancy between the dust temperatures estimated by 
AKARI and those estimated by including submillimetre data. The 
dust mass in M 81 deri ved by the Herschel observation (3.4 x lO'^ 
MM; lBendo et al.ll2010l) is similar to that estimat ed by the AKARI 
obser vation adopted in this paper (3.2 x 10 Mq ; ISun & Hirashital 
l201lh . 
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Figure 9. Same as Fig.[8]but for /in = 0.01. 
5 DISCUSSION 

5.1 Effects of grain size distribution 

We have shown that the grain size distribution significantly af- 
fects the evolution of dust mass above a certain metallicity level 
where the grain growth in clouds is activated. In fact, as shown in 
Figs. [8] and [9l the difference in the grain size distribution makes 
a significant imprint in the relation between dust-to-gas ratio and 
metallicity. This comes from the dependence of /? on metallicity. 
As shown in equation J31t . /3 oc (rgrow/Tci)^^. Fig.|6]shows that 
Tgrow/Tci is a decreasing function of metallicity. Thus, /? increases 
as the system is enriched with metals, and if the last term in equa- 
tion l |48l ) becomes positive at a certain metallicity, the dust mass 
grows in a nonlinear way. Because /? (or Tgrow/Tci) is sensitive to 
the grain size distribution, the resulting relation between dust-to- 
gas ratio and metallicity depends largely on the grin size distribu- 
tion. In other words, the metallicity level at which the grain growth 
in clouds is activated varies sensitively with the grain size distribu- 
tion. This metallicity level is called 'critical metallicity', which is 
quantified and discussed in the next subsection. 

5.2 Critical metallicity for the grain growth 

lAsano et al. I ( l201lh and llnoud boilh have shown that the grain 
growth by accretion dominates the grain abundance if the metal- 
licity is larger than a certain 'critical metallicity'. As mentioned in 
the previous subsection, the nonlinear increase of dust-to-gas ratio 
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Table 3. Critical metallicities for grain growth. 
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^ Grain size distribution: "p" for the power law and 
"(5" for the 5 function. 



by grain growth is realized when — /3sn — TZ + P/e > in equation 
( 148b . Since /3sn ^ TZ, the critical metallicity, Zcr, can be estimated 
by the metallicity at which /3 as a function of metallicity realizes 

/3(Zcr) = e/3sN, (49) 

where P{Z) indicates that /3 is a function of metallicity. In our mod- 
els, € = 0.1 and /3sN = 9.65, i.e., /3(Zcr) = 0.965. Therefore, if 
the metallicity is so high that the grain growth in individual clouds 
increases the dust mass by 0.965 times the original dust amount, the 
grain growth becomes prominent. The critical metallicity varies by 
a factor of ~ 1.5 if we change /3s n by a factor of 2 (i.e . uncertainty 
in /3sN from the assumed species; Ijones &Nuflill201lh . 

The position for the critical metallicity is marked on each line 
in Figs. [8] and Hand listed in Table[3](for /in = 0.1 but Zcr only 
decrease by ~ 10-20 per cent for /in = 0.01). Indeed the dust- 
to-gas ratio rapidly increases if the m etallici t y beco mes larger than 
the critical metallicity. This supports llnouel OOllh's view that the 
activation of grain growth is driven by the metallicity. In this paper 
we have found that the critical metallicity is sensitive to the grain 
size distribution. If the major part of the grains are as large as ~ 
0.1 \x.m, the grain growth is activated above ~ 2 Zq, which is 
too large to explain the large grain abundance in the objects with 
sub-solar metallicities. If the typical grain size is as small as ~ 
0.01 \un or the grain size distribution is described by a power law 
with r > 3.5, the large grain abundance in sub-solar metallicity 
galaxies can be naturally explained by the efficient grain growth in 
clouds. Therefore, the evolutionary history of grain size distribution 
is important to understand the grain abundance in galaxies. 

By an analysi s of infrared dust emission spectra, 
iGalliano et al] ( l2005h show that the grain size distributions 
of some metal-poor dwarf galaxies are biased toward smaller 
grains compared with the Galactic case (i.e. r ~ 3.5). Therefore, 
the critical metallicity for these dwarf galaxies is expected to be 
lower than their metallicities: in other words, the grain growth in 
clouds is expected to contribute significantly to the dust abundance 
in these galaxies. It is interesting to point out that r — 4.5 explains 
the data points of metal-poor galaxies better than r = 3.5. This is 
consistent with Galliano e t al.l iZQOj's conclusion that the grain 
size distribution is biased to small sizes. 

After the grain growth, all the lines in Figs. |8] and |9] finally 
converge. This is explained as follows. If the metallicity is high 
enough, the grain growth is regulated by the lifetime of clouds 
and is independent of grain size distribution. In other words, /3 ~ 
^/(l — ^) if the metallicity becomes high enough (i.e. y is large 
enough in equation[36ll. Using ^ — 1~ D/Z and equation l |48t , we 
obtain the following estimate in the case where the grain growth 
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dominates the increase of the dust content: 



/3s N 



V 

1 i 



v_ 



(50) 



which is vahd for Z > (we only adopt the dominant terms). 
Therefore, if ^ is so small (large) that the right-hand side of this 
equation is positive (negative), ^ tends to increase (decrease) as Z 
increases. This means that ^ tends to approach the value that makes 
the right-hand side to be zero; that is. 



for Z » Zc: 



(51) 



1 + e/?SN 

For the values adopted in our models (e = 0.1 and /3sn = 9.65), 
^ ~ 0.49. This means that the fraction of metals in dust phase is 
about 0.5 at metallicities much higher than the critical metallicity. 
If we consider the uncertainty in /3s n b y a factor of 2 because of 
material difference jjones & Nut"hll201 ih . ^ ~ 0.33-0.66. 

It is interesting that ^ at Z S> Zcr only depends on e/3sN. 
Using the definition of /3sn = e^M^'y /tp (below equation l44l( and 
equation ( |37| l. we obtain 



e/3sN 



(52) 



This is the ratio between the gas mass swept by SN shocks per unit 
time multiplied by the efficiency of dust destruction in SN shocks, 
and the foiTnation rate of clouds. In other words, this is the ratio be- 
tween the dust destruction rate by SN shocks and the grain growth 
rate in clouds. Therefore, it is natural that the final fraction of met- 
als in dust phase can be described by the balance between the dust 
dest ruction by SN shocks and the dust formation in clouds. 

Ilnoud mi l[) also obtained a similar expression for the crit- 
ical metallicity. We have confirmed his conclusion that the dust- 
to-metal ratio app roaches a con stant value: this behaviour is called 
self -regulation by Ilnoud ( 1201 ll) . The difference between our for- 
mulation and his is that he treated rsp as a parameter independent 
of Tgrow while we connect these two parameters through the star 
formation efficiency e. Since both star formation and grain growth 
occur in molecular clouds, these two processes are not independent. 
Moreover, if the grain growth is efficient enough, the dust growth 
time-scale is limited by the lifetime of clouds, which is indepen- 
dent of the grain size distribution. Therefore, the dust-to-metal ratio 
does not depend on the grain size distribution for Z 3> Zcr, and it 
depends only on e/3sN. 



5.3 Significance in galaxy evolution 

The first dust should be produced by SNe in the death of massive 
stars. Because of shock destruction in SNe, the dust sizes may be 
biased to a > 0.1 [x m (Bianc hi & SchneideJl2007l : iNozawa et al.l 
[2007). Hi rashita et a l. (2010) show that shattering driven by inter- 
stellar turbulence can produce small grains efficiently if the metal- 
licity becomes higher than a critical value (~ 0.1-1 Zq). Thus, 
the critical metallicity for the grain growth is near the metallicity 
level where shattering produces small grains efficiently. The pro- 
duction of small grains accelerates the grain growth by accretion, 
which raises the grain abundance. With the increased dust abun- 
dance, shattering can be further efficient, although the produced 
small grains by shatte ring does not nece ssarily activate further shat- 
tering of large grains djones et al .Il996l) . Such an interplay between 



shattering and accretion may be interesting to investigate in the fu- 
ture. 

The size distribution of grains produced by AGB stars is 
also important since AGB stars become the dominant dust pro- 
duction source after several hundred s of Myr dValiante et al .11200^ : 
lOall et aill201 lal : lAsano et al.ll201 ih . The size of grains produced 
in AGB stars is suggested to be large (a ~ 0.1 ^.m) from the 
observations of spectral energy distributions (Groene wegen|[l997l : 
Gange r et al, 1999) , although .Hofmann et al. (2001) show that the 
grains are not single-sized. To clarify the grain size distribution 
foiTned by AGB stars is important for the efficiency of grain 
growth. Shattering in the ISM may also play a role in efficiently 
producing small grains even if AGB stars only produce large grains 
(lHirashitall201Q) . In this case, efficient grain growth can occur. 

It is interesting to point out that the critical metallicity is 
within the metallicity range typical of dwarf galaxies (Table [3](. 
This confirms the conclusion by Asano et al. (201 1) that the strong 
metallicity dependence of the dust-to-gas ratio in dwarf galaxies 
can be explained by the grain growth in clouds. 

The grain growth is also important in some high-redshift 
galaxy populations. In fact, high-redshift quasars have solar (or 
more) metallicities (e.g. Juarez et al. 2009), which implies that the 
grain gr owth is indeed governing the dust abundance in distant 



grain gr owth is maeed governing tne dust abunaance in gistant 
quasa rs dMichalowski et al~ 2010l: [pTpino et ah 201 ll : lAsano et al.l 



20111) (but see lVaUante et alJl200a : lGall et al.ll2011bl) . If the abun 

dance of small grains are enhanced because of shattering as sug- 
gested by Hirashita et al. (2010), the critical metallicity becomes 
lower, so that the importance of grain growth in clouds is further 
pronounced. 

Qualitatively it may be predicted that the mid-infrared emis- 
sion from very small grains is relatively suppressed if the grain 
growth in clouds is activated. However, it is hard to quantitatively 
predict the galaxy-scale observational features caused by the grain 
growth in clouds because it is difficult to selectively see the clouds, 
where grain growth is occurring. Since observations of galactic 
spectral energy distribution inevitably include the emission from 
diffuse medium, other mechanisms modifying the grain size distri- 
bution such as shattering and coagulat ion are also reflected in the 
observed emission from grains. Indeed, lGamanoeta"n ( l2005 l) show 
that the grain size distribution is biased toward smaller grains in 
some dwarf galaxies, which may be interpreted as the efficient grain 
processing in diffuse medium. They also demonstrate that there is a 
variety in the grain size distribution among dwarf galaxies. Such a 
variety will be investigated in the future with a consistent treatment 
of a nonlinear combination between the small grain production by 
shattering and sputtering and the grain growth by accretion and co- 
agulation. 



6 CONCLUSION 

We have formulated and investigated the grain growth rate by ac- 
cretion in interstellar clouds. The formalism is applicable to any 
grain size distribution. We have found that the grain size distribu- 
tion is really fundamental in regulating the grain growth rate. We 
have also implemented the formulation of grain growth in individ- 
ual clouds into the chemical evolution models of entire galaxies. 
The models also treat dust supply from stellar sources and dust de- 
struction by SN shocks, but we have focused particularly on the 
grain growth in clouds in this paper. We have found that the metal- 
licity level where the grain growth in clouds becomes dominant 
strongly depends on the grain size distribution. If the significant 
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fraction of the grains have radii < 0.01 (xm or the grain size dis- 
tribution is described as power law with r > 3.5, the large grain 
abundance at the sub-solar metallicity level is naturally explained 
by the grain growth in clouds because the surface-to-volume ratio 
of the grains is large enough. The grain growth should be efficient 
in galaxies whose metallicity is above Zcv estimated in Section [5!2l 
Our formulation for the grain growth is applicable to any grain size 
distribution and is implemented straightforwardly into any frame- 
work of chemical enrichment models. 
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APPENDIX A: INSTANTANEOUS RECYCLING 
APPROXIMATION 

We estimate TZ and yz used in Section l44l With an initial mass 
function (IMF) (j>{m),TZ and yz are written as 



TZ = I [m — w{m)]<l){m) dm, 

J rat 

yz = mpz{m)(l>{m)dm, 



(Al) 
(A2) 



wiiere rrit is tlie turn-off stellar mass, rriu is the upper mass cutoff of 
stellar mass, m is the stellar mass, Wm is the remnant mass, pz {m) 
is the fraction of mass converted into metals in a star of mass m. 
We assume the Salpeter IMF ((^(m) cx m~^ "^^) with stellar mass 
range 0.1 Mq ^ m ^ 100 Mq. The IMF is normalized as 



m<j){m) dm = 1, 



(A3) 



where m\ is the lower mass cutoff of stellar mass. 

For t he rem nant mass, we adopt the fitting formula provided 
bv llnoudJioTTI) : 



w(m) 



= < 



r 1 

0.13 
0.13 



8M(7 







(m > 40 Mq), 
(8 ^ m ^ 40 M(7 



(A4) 



8M(7 







(m < 8Mq). 



We also adopt the fitting formula for the mass of ejected metals as 
a function of stellar mass as tlnoueii2011.) 



mzvm) 



m 



= < 





0.02 
0.02 



8Mq 



8Mq 



(m > 40 Mq), 
(8 ^ m s: 40 Mq) 

(m < 8Mq). 



(A5) 



llnoud dm l|) show that there is no trend with metallicity Z for 
mz/m. In fact, raz is the mass of ejected metals, so we need to 
subtract the metal mass already included before the nucleosynthe- 
sis: mpz(m) = mzim) — Z[m — ui(m)]. Since the metallicity 
range of interest is ^ ^ < 0.02, mpz{m) is reasonably between 
mz{m) — 0.02[m — w(m)] (called minimum) and mz{m) (max- 
imum). In Fig. lAll we show TZ and yz as a function of the turn-off 
mass (or age). We adopt the values at t = 5 Gyr, that is, TZ — 0.18 
and yz = 0.013 (average of the maximum and the minimum), 
since the typical gas consumption time-scale or sta r formation time- 
scale of nearby star -forming galaxies is 1-10 Gyr dKennicutll 1 9981 : 
iGavazzi et alj2002l) . However, the uncertainty caused by the age is 
within a factor of 2 if we adopt an age of > 4 x 10^ yr. Thus, as 
long as we treat galaxies whose typical star formation time- scale 
is longer than a few x lO'^ yr, our calculation gives a reasonable 
results for the metal and dust enrichment. 

The SN rate, 7, is also necessary to estimate the dust destruc- 
tion rate. If we assume that stars with m ^ 8 Mq become SNe, 
7 can be related to the star formation rate by approximating the 
lifetimes of these stars to be zero: 
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Figure Al. The returned fraction of gas TZ (solid line) and the fraction 
of newly produced metals yz (dotted and dashed lines) are shown as a 
function of turn-off stellar mass or age. For yz, the maximum (dotted line) 
and the minimum (dashed line) cases are shown. 
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